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Abstract 



> 

^ ■ The Wheeler-DeWitt equation is investigated and used to examine a state 

o 

Q\ • after a quantum tunneling with gravity. To make arguments definite we treat 

o. 

a discretized version of the Wheeler-DeWitt equation and adopt the WKB 
method. We expand an Euclidean wave function around an instanton, by us- 

O" 1 

5h . ing a deviation equation of a vector field tangent to a congruence of instantons. 

The instanton around which we expand the wave function corresponds to a 
so-called most probable escape path (MPEP). It is shown that, when the wave 
function is analytically continued, the corresponding state of physical pertur- 
bations is equivalent to the vacuum state determined by positive-frequency 
mode functions which satisfy appropriate boundary conditions. Thus a quan- 
tum field theory is effective to investigate a state after a quantum tunneling 
with gravity. The effective Lagrangian describing the field theory is obtained 
by simply reducing the original Lagrangian to a subspace spanned by the 
physical perturbations. The result of this paper does not depend on the op- 
erator ordering and can be applied to all physical perturbations, including 
gravitational perturbations, around a general MPEP. 
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I. INTRODUCTION 



The universe is thought to have been experienced phase transitions many times. Some of 
these phase transitions may be due to quantum tunneling. Hence there arises an intellectual 
and practical urge to know a state after a quantum tunneling. For this purpose many authors 
investigated the Schrodinger equation of scalar fields in flat spacetime |T]J^] in the multi- 
dimensional tunneling approach ||. They showed that a state after a quantum tunneling in 
flat spacetime is the vacuum state determined by positive-frequency mode functions which 
satisfy appropriate boundary conditions. The result says that a quantum field theory of 
the scalar fields is effective to seek a particle spectrum after a quantum tunneling if the 
boundary conditions are attached on the mode functions. The effectiveness makes it possible 
to calculate a CMB anisotropy of the universe theoretically in some inflationary scenarios [|J . 

However the calculated anisotropy is based on the formalism developed in flat spacetime, 
as mentioned above. Since gravitational effects can be significant in some cases , we want 
to generalize the formalism to include gravity. As far as only the scalar fields' perturbations 
are concerned, the generalization was done in Ref. ||, by neglecting gravitational pertur- 
bations. In this case, a state of the scalar fields' perturbations after a quantum tunneling 
with gravity is obtained by the same method as in flat spacetime, except evolution equations 
of mode functions is replaced by those on the curved background. Thus a quantum field 
theory of the scalar fields in curved spacetime is effective to investigate a state of the scalar 
fields' perturbations after a quantum tunneling, provided that gravitational perturbations 
are neglected. However, when we intend to investigate a state of gravitational perturbations 
after a quantum tunneling (for example, when we calculate a spectrum of the primordial 
gravitational waves), we have to include gravitational perturbations. The corresponding 
formalism must be derived from a more fundamental level. 

When we investigate a classical dynamics of gravitational perturbations, so-called gauge 
invariant variables play a important role. The gauge invariant variables are those linear 
combinations of perturbations which are invariant to linear order under coordinate trans- 
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formations, and a maximal set of them describes all physical degrees of freedom of linear 
perturbations. For several background geometries, it was shown that classical dynamics 
of the gauge invariant variables can be described by some auxiliary scalar fields (Ref. |7j 
for Minkowski background, Ref. [|8] for the Robertson- Walker background, Ref. for the 
Milne background, and Ref. []TU[ for the one-bubble inflationary background). Hence we can 
expect that a quantum field theory of the auxiliary scalar fields can be applied to investigate 
a state of gravitational perturbations after a quantum tunneling, although the equivalence 
is established only in a classical level. 

At this point, we can expect that a quantum field theory will be effective to investigate a 
state of all physical perturbations after a quantum tunneling with gravity. Thus we want to 
show from a quantum theory of gravity that the expectation is true. In all papers with this 



aim 



nfl ) the Wheeler-DeWitt equation is expanded by some collective coordinates around 
a background. Our strategy in this paper is to expand not the Wheeler-DeWitt equation 
but a wave function itself covariantly. Anyway, due to the result of this paper, we can safely 
use the quantum field theory of perturbations in order to calculate a CMB anisotropy, a 
spectrum of the primordial gravitational waves, etc. in some inflationary scenarios. 

Now let us review a quantum theory of gravity based on a canonical formalism. Con- 
sider a system of a n-dimensional spacetime (M, g^ v ) without boundary and scalar fields \ a 
described by the action 

1 „ 1 



Ie-kg = J d n x^g 



-- 2 R - ^l ab (x c )d, X a d„X b - V( X ° 
2k z 2 



where ^ab is an arbitrary positive-definite matrix. In the ADM decomposition |7j of the 
spacetime 

ds 2 = -K\dt 2 + q ij (dx i + A i dt)(dx j + A j dt) , (1.2) 

functional derivatives of the action Ie-kg with respect to Aj_ and A 1 give the following 
constraints. 

H±x=0 , 
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ft ia; = -2^ + 7r a 9 iX a 
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The former is called Hamiltonian constraints and the later is called momentum constraints. 
In the expressions, n 1 ^ is a momentum conjugate to g^, n a is one conjugate to x a , Dj denotes 
a covariant derivative compatible with g^-, and 



for 2, 



„ 1 / 2 

(j-jjfei = — — yqiklji + <?ii<?jAi — — — ^QijQkl 



Gijki = —= (qikQjl + QilQjk) 
v ^ 



;i.5) 



;i.e) 



for n = 2. It is well-known that Poisson brackets among the above constraints are given by 

{fH ± , gU L } P = (fD l g - gD l f)H t , 
{f i H i ,gH ± } P = (f%g)n ± , 

{fHi, g J 'Hj\r = (pl),g' - g-'I^fYH, , (1.7) 

where fH± and f z Td denote / dxf(x)Tl±x and / dxf t (x)Hix, respectively. These show that 
the constraints are first class. A method to quantize a system with first class constraints was 



given by Dirac |12| and leads the following simultaneous differential equations for a wave 
functional 



H± X V = 



;i.8) 
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where H_lx and Tiix are differential operators obtained by replacing 7T lJ with —ih5/5qij, n a 
with —ih5/5x a in H±x and 7iia;, respectively. Equation fll.Sfl is called the Wheeler-DeWitt 
equation |13| . In defining the differential operators, there is a problem of operator ordering. 



The first principle to determine the operator ordering is requiring the following algebra [14 



Wu 0H±\ ^ = iHfD'g - gD'fiHi* , (1.9) 

y%,gH±\y = iHf l dig)H±y , (l.io) 

f i n i , g m j ]^ = ih[f,g} i n^ (1.11) 

for an arbitrary functional ty. The algebra is a quantum version of ( |1.7| ), and is called the 
Dirac algebra [|l2j. The condition (|1.9Hl.llj ) is necessary in order for the equation (|1.8|) to 



give a consistent quantum theory. 

The remaining part of this paper is organized as follows. In Sec. || we investigate a 
simultaneous differential equations, which become the Wheeler-DeWitt equation in a limit. 
We expand an Euclidean wave function around an instanton, by using a deviation equation 
of a vector field tangent to a congruence of instantons. The instanton around which we 
expand the wave function corresponds to a so-called most probable escape path (MPEP). 
In Sec. |HT| , by analytically continuing the wave function, it is shown that the corresponding 
state of physical perturbations after a quantum tunneling is the vacuum state determined 
by positive-frequency mode functions which satisfy appropriate boundary conditions. Hence 
a quantum field theory is effective to investigate a state after a quantum tunneling. Sec. [TV] 
is devoted to summarize this paper. 



II. EUCLIDEAN WAVE FUNCTION FOR A DISCRETIZED 
WHEELER-DEWITT EQUATION 

Let {M.xi Gxap) be a family of D-dimensional pseudo-Riemannian manifolds parame- 
terized by an integer x (= 1, 2, ■ ■ • , x max ). A configuration space we consider is a pseudo- 
Riemannian manifold constructed from them as 



5 



(M, G a p) = ® X {M X , G xa p) . (2.1) 

Then we consider the following simultaneous differential equations for a complex function 
^ on M. 

H± x V = , 

H I ^ = , (2.2) 

where x = 1, 2, • • ■ , x ma:r , J = 1, 2, ■ • • , J maa; and 

ft 2 

H± x = — —GlfD a Dp + V x , 
W/ = -ihv?D a (2.3) 

are linear differential operators on Ai. As usual we call a solution of this equation a wave 
function. In the expression D is a covariant derivative compatible with the metric G a p, 
{V x } is a set of functions on Ai, and {f?} is a set of linearly independent vectors on 
A4. The Wheeler-DeWitt equation ( |1.8| ) can be written in this form provided that it is 
properly discretized and the operator ordering problem is ignored. Ti,± x corresponds to the 
Hamiltonian constraint at the point x on a spacelike hypersurface, and Tii corresponds to 
the momentum constraint for a point and a direction both of which are specified by /. For 
example, 

Q<hj{x')q k l{x") = 2 K 2 G ljkl (x)5 xx ,5 xx „ , 
Qqij{x') X a (x") = q 



gtfVjxV') = ^^^^^ {2A) 
'q(x) 



for the system considered in Sec. |, where Gijki is defined by ( |1.6| ). We mention that terms 
linear in the derivative, which may appear when the operator ordering problem is solved, 
were not included in 7i± x for concreteness, since they do not change our conclusion. 

For the equations to be solved consistently it must be assumed that any commutators 
between the linear operators can be written as linear combinations of themselves when they 
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Hi,Hj\ ^ = ihcf) K H K y 



operate on an arbitrary function on the configuration space. Equivalently they must generate 
a Lie algebra with respect to the commutators. In this paper we adopt the following algebra 
of commutators between the differential operators. 

H ±x ,H^]y = thc^H I y , (2.5) 
Hi,H± x ]* = ihc?] x ' Wiv* , (2.6) 

(2.7) 

for an arbitrary function \1/ on M.. Note that the right hand side of ( |1.9|) , ( 1.10|) and 
( |1 . 1 1| ) are linear in 7i!j\I/'s, 7ij_\l/'s and 7Tlj\I/'s, respectively while the right hand side of ( |2.5| ), 
( |2.6| ) and Q2.7| ) are linear in 7i/\l/'s, H^^'s and 7l!/\I/'s, respectively. Hence the algebra 
( p.5| - [2.7p is a generalization of the discretized version of the Dirac algebra ( |1.9[ - |1.11] ). We 
include the case when the 'structure constants' and c^ 3 ' depend on a position in 

the configuration space M.. Although confirmation of the algebra requires a knowledge of 
short distance behaviors of the theory, we simply assume that the algebra does hold. After 
detailed investigation of the discretized system, we will extract results independent of not 
only the operator ordering but also a way of the discretization. Note that since a dicretized 
version of the Dirac algebra (|1.9| - |l-llD can be written in the form (|2.5| - |2.7|) , the system we 
consider includes the discretized Wheeler-DeWitt equation as an important example. Now, 
for later conveniences we rewrite the assumption ( |2.5| - 12.7| ) as follows. It is equivalent to the 
following set of equalities 

GfdpV x ,-G Q fdpV x = -c { ^vf , (2.8) 
G?D n 4 + G*Drf = cfl x 'G a J , (2.9) 
v?d a V x = -c { j% x 'v x , , (2.10) 



[v I} vj] a = -c 



IJ V K 



(2.11) 



and 



GfD a D p V x , - G a jD a D p V x = , 

GfD a D p v] = -GfvfR\ sp , 



(2.12) 
(2.13) 



where R is a curvature tensor of D: 



R{X,Y)Z = D X D Y Z - D Y D X Z - D [xy] Z . (2.14) 



Note that, while ( |2.12| - |2TT3| ) are affected by the operator ordering, ( |2.8| - |2.11| ) are not. In 



particular the fact that {vf} is integrable is independent of the operator ordering since it is 
a consequence of Q2.11|) . We denote the integral surface G. In this paper the quotient space 



Ai/G plays a important role since the second equation of (|2.2| ) shows that G has nothing to 
do with physical degrees of freedom. In the case of the usual Wheeler-DeWitt equation the 
quotient space A4/G is called a superspace. Hence we call the quotient space a super space in 
our case, too. Evidently, the superspace is a space of all physically distinct configurations. 

Next we investigate a classical mechanical system whose quantum version corresponds 
to the equations ( |2.2|) . Let us consider a phase space (T*Ai, dp a A dq a ) and a Hamiltonian 
of the form [] 

H = A^H± X + A T Hi , (2.15) 
where (7^ 0) and A 7 are Lagrange's multipliers and 

n ±x = ^Gf PaP/3 + v x , 

Ui = v a lVa . (2.16) 



In the expression, p a (e T*M) is a momentum conjugate to a coordinate q a of M.. (|2.8[ - |2.1lD 
is equivalent to the following algebra of Poisson brackets. 

{H± X} H± X '} P = c^fHi , 
{Hi,H± x } p = cf^ x H± x > , 
{Hi,Hj} p =cfj K H K ■ (2.17) 



1 Hereafter we apply the Einstein's summation rule unless otherwise stated. 



S 



The equivalence is consistent with the well-known fact that the operator ordering does not 
affect the corresponding classical dynamics. The classical equations of motion in the config- 
uration space M. is obtained from the constraint equations and the Hamilton's equation as 
follows 0. 

-G^JV^ + (Al) 2 V X = , 

J\f l3 T>pN' a = —Q^dpV + Q a ^v J (dpWy — d^wp) , (2.18) 



where 



J\f a = ^ - v a , 
dt 

v a = A*vf , 



g af3 = A x ± Gf , 



Qafi = \ Q 



>-l 



a/3 ' 



V = A X ± V X , (2.19) 



and 



w a = g a ^ . (2.20) 

In the expression T> is a covariant derivative compatible with the metric Q a p. 

Now we return to the problem of solving the simultaneous differential equations ( p.2| ) to 
obtain a wave function. To extract those properties of the equations which are independent 
of the operator ordering we adopt the WKB method. First, without loss of generality, we 
can expand the wave function \1/ as 



^ = exp 



1(W {0) + HW {1) + ■■■) 



(2.21) 



where W {0 \ W {1 \ ■■ ■ are complex functions on M.. Next we solve the differential equations 
order by order in h, considering h as a small parameter. 



Here the Einstein's summation rule is not applied with respect to x. 
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A. Lowest-order WKB wave function 



To the lowest order in h the differential equation ( |2.2| ) is reduced to the following simul- 
taneous differential equations for . 

\cf d a W m dpW {a) = V x , 

v?d a W {0) = . (2.22) 

These determines as a function on A4, provided that a suitable boundary condition is 
attached. When W^ ' is real in a region of the configuration space, we call the corresponding 
wave function \1/ an Euclidean wave function. We call the region an Euclidean region. In 
the remaining of this paper we investigate the Euclidean wave function. Physical meaning 
of the Euclidean wave function will become clear in the following arguments. 

For the real function we can define a family of vector fields {N£} parameterized by 
the integer x: 

N« = GfdpW {0) . (2.23) 



Operating the differential operator G^fDp on the first equation of ( |2.22|) , we obtain the 
following equation. 

N x D P N x> = G°fdpV x , (2.24) 

where we have used the fact that the covariant derivative D is torsion free and compatible 
with G xa p. Note that the compatibility with G xa p is a result of the direct-sum structure 
(0) of (M,G a p). 

Since the first of ( [2.22 ) is of the form of the Hamilton- Jacob i equation, we may regard 



Q2.24Q as the corresponding equations of motion for the system. In fact, if we define a vector 
field Af a by 

M a = A^iV" , (2.25) 

then we can show that 
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N p V p N a = Q ap dpV . (2.26) 

Since dpu y — d^up = in this case, these equations show that Af a satisfies the classical 
equations of motion (|2.18 ) in which V x is replaced by —V x , where u a = Q a pN . Therefore 



M a defined by ( |2.25|) is a tangent vector field of a congruence of instantons. In this sense 
the Euclidean wave function define a congruence of instantons. Conversely, if a congru- 
ence of instantons satisfies the condition v^^dpw^ — d^wp) = 0, then we can construct the 
corresponding lowest-order Euclidean wave function as follows. Introduce a parameter r by 

and 

= 2V . (2.28) 



B. Expansion of the Euclidean wave function around an instanton 



In this subsection we consider a real solution of ( |2.22|) and expand it around an instanton. 



An expansion of the corresponding Euclidean wave function is obtained from that. 

The second of the constraint equations ( |2.2|) or the second of ( [2.22|) suggests that any 



equations for physical degrees of freedom can be written in forms invariant under any dif- 
feomorphism of G. Hence we can expect that there is a set of real functions {A^} such that 
a weighted metric X]_G^ and a weighted 'potential' \±V X are invariant under the diffeo- 



morphism. We first investigate the weighted metric. With the help of ( |2.9| ) it can be shown 
that the condition L„ / (A^G° /3 ) = is equivalent to 



v ( }d a Xl = Xic { ^ , (2.29) 
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where L represents a Lie derivative in the configuration space A4. The invariance of the 
weighted 'potential' L„ / (A^14) = is equivalent to ( |2.29| ), too. Therefore both of \ x ± G°f 
andA^\4 are invariant under the diffeomorphism if and only if {A^} satisfies ( |2.29| ). What 



is the meaning of the condition ( [2.29P ? cf]f is the 'structure constant' which appears in 
the algebra ( |2.5| - |2.7D and represents a way of transformation of the constraint 1H± X under a 
group generated by {Hi}, which corresponds to G. Hence the condition (|2.29| ) means that 



the coefficient A^ must change covariantly under the group transformation. In order for 
the equation (|2.29|) to be solved consistently, the following integrability condition must be 



satisfied: 

(3)K (2)x , (2)x" (2)x _ (2)x" (2)x , a o (2)x _ af) (2)x _ Q /„ 
C 7"J Kx 1 Ix' c Jx" c Jx' c Ix" < u I u a c Jx' u J u a c Ix' ~ u • ^.OUJ 



It can be easily confirmed that ( |2.3U| ) is actually satisfied as a consequence of ( 2.10|) and 



( p. 11| ) 0. Let {A ± ^/} be such a complete set of linearly independent real solutions of (|2.29|) 



that 

^±x'\£I=£l = &x' > (2-31) 

where {£ 7 } is a coordinate system of G and {£q} is a set of constants. Then a general 
positive solution of (|2.29j) can be written as a linear combination of these solutions: 



A1=AIA ± :, , (2.32) 

where {A^} is a set of positive functions on the superspace M./G in the sense that it is a 
set of positive functions on M. satisfying \j Vl h 3 \_ = 0. In the remaining of this subsection we 
restrict A^ to this form. 

From the set of coefficients {Aj_^/} we can define a set of 'metric' {G"^ 13 } and a set of 
'potentials' {V x } by 



3 The consistency condition (2.30) can be understood as a consequence of Jacobi identities derived 
from the commutators (fHj-2,7). 
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^ x — A ±x ^x' 1 

V X = \JV X , . (2.33) 

Here both of V x and G'% are invariant under the group transformation of G. However, since 
G"^ 3 may have components in the direction of G, it can not be regarded as a tensor on the 
superspace M.jG without any modification, while V x can be. We want to modify it in order 
to regard it as a tensor field on M./G. For this purpose define the following 'projection 
operators'. 

nad r" a a IJ 

gap _ g y _ v uyj v r ^ 

Ga0 = G a^G^G u0 ) 

y p = y y u0 , 

gf = g' a ^ , (2.34) 



where 



G' a p = [G'- 1 ) a0 , 

r>i a 

ll3 = G a0V I V f } , 

13- (^ IJ 



1 U = (l- L ) ■ (2.35) 
It can be easily confirmed by using (2.11) that 



h Vl Q aP = h Vl G a p = L Vl G a = KA P = • (2-36) 

Then define the modified tensor G^ by 

Gf = G a ,G'^G/ . (2.37) 

It is evident that G^ can be regarded as a tensor field on the superspace A4/G since 
L Vl G x P = and G^Gp-yv] = 0. Note that G alS equals to a weighted sum of G^ as 

ga0 = XlGf , (2.38) 
13 



and it can also be regarded as a metric tensor of the superspace M./G. In terms of the 
metric tensor we can rewrite the ( |2.22| ) as the following simultaneous differential equations 
on Ai/G for a real function on Ai/G. 

-Gf daW^dpW^ = % . (2.39) 



Our final goal in this subsection is to expand a solution of (|2.39|) around an mstanton 



in the superspace. For this purpose we consider such a congruence of instantons in the 
superspace that each instanton is an integral curve of the vector field M a defined by 

M a = g al3 d^W (0) , (2.40) 

and introduce an Euclidean time f as a coordinate variable in Ai/G by 

• (2.4D 

If we define a set of vector fields {N£} by 

NZ = Gfd W<® , (2.42) 

then 

M a = AlN« , (2.43) 

and N" generates a time reparameterization. Commutators between N£ and N£ is zero 
by using fl2.8|) and the fact that the covariant derivative D is compatible with G°f . Hence 
the set of vector fields {N£} is integrable in A4/G. Then a coordinate system {f, rj x , ip n } 
of Ai/G can be introduced so that {f,r] x } is a coordinate system of the integral surface 
of {N£}- For these definitions, dependence of on {f,r] x } is completely determined by 
Q ) as 

N«d a W {0) = 2V X . (2.44) 

Note that there is a one-to-one map from {tp n } to a space of all integral surface of the set 
of vector fields {N®}- Moreover {tp n } is a maximum set of coordinate variables on which 
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W^°)'s dependence is not completely determined by (|2.39|) . Thus we can say that the set 
of coordinates {(p n } represents all physical degrees of freedom of perturbations around the 
instanton. The coordinate f represents an Euclidean time of the congruence of instantons, 
and {r] x } represents degrees of freedom of time reparameterizations. 

Now let us investigate the dependence of on the physical coordinates {y? n }. For 
this purpose we derive a deviation equation of a vector field in Appendix [A|. Applying the 
resulting deviation equation to the vector fields M a and = (d/dip n ) a in the superspace 
(Ai/G, Gap), we obtain 

V F1 Z% = AQ^VpW® , (2.45) 

v FL z% = L zi L g a ~< (v^VpV - % pfia MpM° - AayVfyv) , (2.46) 

where T> is a covariant derivative compatible with the metric G a p, f)p denotes a Fermi 
derivative along Af a made from T>, and 

= Zl- (Q(Af, Zn)/G{N,N)) M a , 

H{X,Y)Z = V X V Y Z -V Y V X Z -V [x ,y]Z . (2.47) 



The second equation ( |2.46| ) makes it possible to investigate how the vector field ±Z£ evolves 



along a integral line of M a \ it corresponds to a linearized equation of motion for the physical 
degrees of freedom of perturbations around an instanton. With the help of the evolution 
equation of perturbations, then, the first equation ( |2.45| ) says that the second derivative of 



in the direction of {±Z?} is determined as follows. 

± Z^ ± ZiV a VpW^ = ± Q a p ± Z«V F± Zi . (2.48) 

Using these results we can expand a real solution of the guide equation around an 
instanton. For this purpose let us introduce a set of new variables {±tp n }, each of which is 
an affine length of a geodesic normal to the instanton such that ±ip n = along F and 

= ± Z« along T . (2.49) 



d 
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Here T denotes the instanton around which we intend to expand a real solution W^ ' of the 
guide equation. Denote the value of rj x along the instanton r by t]q. Then we can expand 
by ±Lf n around the instanton T as follows. 



[2.50) 



where the matrix Q a p is defined by 



Z- 1 ) V Fi Z2 



(2.51) 



and ±Z£ follows the evolution equation ( 2.46|) . In the expression, (^±Z x ) is defined by 

±Z° ^Z- 1 )* = 6* , (2.52) 

and so that ±2i is a projection operator to a subspace generated by {±Z£}. 

As a result, a general real solution of ( 2.39|) can be expanded as follows. 



W^(f, V l ± r) = W®(f,rf*,0) + ~O a/3 (f) ± r ± / + O ( ± 3 ) 



(2.53) 



where the new variables {±4> a } are defined by 



rya , r. 



(2.54) 



provided that its dependence on {f, rj x } is fixed by ( |2.44| ). Note that the set of variables 
{±<f) a } represents all physical degrees of freedom of perturbations around the instanton. 
From the expansion ( |2.53p we obtain the corresponding expansion of the Euclidean wave 
function. 



#oO")exp 



+ 



(2.55) 



We have expanded a general real solution of ( |2.39| ) around an instanton T. Finally let 
us restrict a class of instantons by imposing a specific boundary condition on the expanded 
wave function. We impose the condition that the matrix Q is positive definite in the subspace 
generated by {±Z£} along the instanton T until a possible turning point, where a turning 
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point means a point on T at which V becomes zero. The positivity of the matrix means 
that probability falls off leaving from the instanton to all physical directions which are 
orthogonal to the instanton. It seems satisfied if the instanton corresponds to a so-called 
most probable escape path (MPEP) [Q. In general, if any first class constraints exist for a 
system, they associate with symmetries of the system and there arises a zero mode problem. 
However, since by definition Q is a matrix on a linear space of only physical degrees of 
freedom, Q has no zero modes which associate with the first class constraints (|2.2|) . Thus 
the positivity is necessary for the instanton to have a meaning of a MPEP. Note that we can 
attach this condition for only a class of instantons, while the expansion Q2.53| ) is possible 
around an arbitrary instanton. That is because the condition on Q may restrict behaviors 
of the corresponding wave function and may violate some physical criterion (eg. asymptotic 
behavior, etc.). Hereafter we restrict a space of all instantons to the class of those for 
which we can attach the positivity of Q, and we adopt the positivity condition. We derive 
significant results from this condition in the next section. 



In this section let us see how a quantum field theory of perturbations is derived from 
the Euclidean wave function expanded in the previous section. First we discuss a continu- 
ation of the wave function through the turning point. After that we give a field-theoretical 
interpretation of the wave function. 



As shown in the previous section, the Euclidean wave function is expanded around an 
instanton T as 



III. QUANTUM FIELD THEORY OF PERTURBATIONS 



A. Analytic continuation of the Euclidean wave function 



= ^o(r)exp 




} 
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where ±<f) a = corresponds to the instanton T. If we regard ±<p a as a quantity of order 
0{Ti p ) (0 < p < 1), then we can neglect the terms of order 0(j_(p a3 , ±4> a h) as a consistent 
approximation: 

1 

2h~ 



*(r,^, ± a ) = vl> o (f)exp 



(3.1) 



It is well known that near a turning point the WKB method is not good. So we need 
a matching condition of WKB wave functions at the turning point. It is one of the most 
difficult problems in the WKB approach to the Wheeler-DeWitt equation. The difficulty is 
mainly due to the following two facts: (1) the configuration space M. is multi-dimensional 
(infinite dimensional in a continuous limit); (2) the 'super metric' Q a p is not positive definite 
in the subspace generated by {±2?}. Vachaspati and Vilenkin |L5| attacked the problem by 
using a simple model and neglecting the second difficulty. They investigated a Schrodinger 
equation in a two dimensional configuration space with a positive definite 'super metric' 
and obtain the following matching condition of lowest-order WKB wave functions for the 



tunneling boundary condition |T6[: if the exponent of the wave function is expanded around 
a classical solution to second order in perturbations orthogonal to it, then the wave function 
of perturbations is continuous at a turning point. Moreover their result insist that a proper 
WKB wave function beyond the turning point is obtained by an analytic continuation. As 
a result the matching problem for their system is equivalent to one for a one-dimensional 
quantum system. Returning to our system, if the weighted 'potential' V does not depend 
on (f n near the turning point with an enough accuracy, then the matching problem for our 
system seems equivalent to one for a one-dimensional quantum system. Hence, in this case, 
we can expect that their procedure does work. Then we obtain the following form of the 



wave function beyond the turning point by the analytic continuation with r — > it. 



$ = W (it) exp 



(3.2) 



where the curve ±4> a = corresponds to that solution of the classical equation of motion 
Q2.18] ) which is the analytic continuation of the instanton T, and {±<f) a } denotes physical 
perturbations around the classical path. Here the matrix Vt a p can be written as 



tt a p(it) = -i±Q aj {±_Z x \v F1 _Z1 , (3.3) 

where T>p denotes a Fermi-derivative along the analytically-continued classical path and 
±2£ is the analytic continuation of ±Z£. 

B. Field-theoretical interpretation 

In quantum field theory, normalized mode functions play a central role. The mode 
functions are normalized in terms of a so-called Klein-Gordon inner product. Thus, for our 
discretized system, we propose the following definition of an inner product between two 
complex vector fields: for complex vector fields X a and Y a on the superspace Ai/G, 

(X, Y) KG = -i ± Q af} ( ± X°"D F± V* - ± Y* a V F± X?) , (3.4) 

where * denotes a complex conjugation. For this definition, it can be easily proved that the 
inner product between {±Z£} is constant along the classical path: 

d - 

77f(_L^m, ±Z n ) KG = . (3.5) 

at 

We want to normalize the set of vectors {±Z£} since each of them is a solution of 
perturbed equation of motion and corresponds to the mode function in a continuous limit. 
By introducing a complex regular matrix whose elements are constant, we take linear 
combinations of complex conjugates of {±Z?} as 

± u a k = Q ± Z*<* . (3.6) 



Evidently is a candidate for the normalized mode function. By using ( |3.3|) we can 
express the inner product between in terms of Q as 

U u k,±Uk')KG = J_U(Q + ^)±U^ ^ . (3.7) 

The positivity of the Q in the Euclidean region, which has been introduced in the last 
paragraph of subsection [QB|, means that the matrix (j_Uk, ±Uk>)KG calculated in ( |3.7|) is 
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positive definite at the turning point. Hence, due to the constancy of the matrix elements 
implied by ( |3.5| ), we can choose the transformation matrix Q so that 

(_|_"Ufc, J_Uk')xG — <W (3-8) 

along the analytically-continued classical path. Then is specified as a set of orthonor- 

malized solutions of the simultaneous equations 

^±«j=-±«f±^^p /9 v+^(^y^ -h 9 ^) ' (3 - 9) 

= Qa/3^ a ±4 , (3.10) 

and is a complete set of basis vectors in a subspace generated by {{d/d±ip n ) a }. Equa- 
tion (|3.9| - |3.10D is a linearized evolution equation of orthogonal perturbations around the 
analytically-continued classical path. The boundary condition of j_u% is that the matrix VL a p 
defined by 

n a(3 = -n_Q ai (±u- l )fv F ^ (3.11) 

is real and positive definite in the subspace generated by {±Z£} when it is analytically- 
continued back to the Euclidean region (see arguments in the last paragraph of subsec- 
tion |I B|) , where * denotes a complex conjugation and (j.^ -1 )^ is defined by 

±<UOa=^' , (3-12) 

and so that (±w _1 )a±wf is a projection operator to the subspace generated by {{d/d±ip n ) a }. 

Finally we show that the wave function ( |3.2| ) can actually be interpreted in terms of 
a quantum-mechanical system, which will become a quantum field theory in a continuous 
limit. Consider a quantum- mechanical system described by the hamiltonian 

H=l(±G aP ±n a ±7c p + ± V al3x } a ± ^) , (3.13) 

and the equation of motion 

ihV F d=[6,H] (3.14) 
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for an arbitrary operator O, where ±if a is a momentum conjugate to ±<p a and 

+ ~^d,Vd u vj . (3.15) 
The equation of motion says that the operator ±<j) a can be expanded as follows. 

J .^ = n 1 / 2 2(a fc±u « + aU<) , (3.ie) 
k 

where a k and annihilation and creation operators which satisfy 



0i k , Q,y 



>kk' 



[a k , a k >] = 



. (3.17) 



For the representation G3.16Q with the normalization ( |3.8| ), the annihilation operator a k can 
be expressed as a linear combination of terms of the form 

±TT a -i^ a0 ±^ . (3.18) 

Hence the wave function ( |3.2| ) is annihilated by all the annihilation operators: 

a k m = (3.19) 



for all k. This says that the wave function (|3.2| ) represents the vacuum state determined 
by the vectors {±u^}. Thus, in a continuous limit, the corresponding wave function can be 
interpreted in terms of a quantum field theory: it represents the vacuum state determined 
by that set of positive-frequency mode functions which is a continuous limit of {±_u%}. 



C. Reduced Lagrangian 

When we intend to apply the formalism developed in this paper, what we have to do 
is: (1) to fix a MPEP; (2) to give a coordinate system in a space of all gauge-invariant 
perturbations (physical perturbations) around the MPEP; (3) to seek explicit form of the 
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effective Hamiltonian (|3.13|) and the matrix (|3.11| ); (4) to calculate any wanted quantities 
by using the quantum field theory described by the effective Hamiltonian. The positive- 
frequency mode function is defined so that the matrix ( |3.11| ) is real and positive definite in 



the subspace spanned by the gauge-invariant perturbations in the Euclidean region. 

In this subsection we show that the effective Hamiltonian and the matrix can be ob- 
tained from an reduced Lagrangean of the gauge-invariant perturbations. Thus the reduced 
Lagrangean gives a convenient prescription to perform the procedure (3). 

The classical Hamiltonian ( p,15| ) corresponds to the Lagrangian 

L = p a q a - H 

= \Q a p(q a -v a )(f-v' 3 )~V . (3.20) 

This Lagrangian describes a dynamics in the full configuration space M.. When we go to 
the superspace M/G, the corresponding reduced Lagrangian L in M. /G is 

L = \Gattnf-V , (3-21) 

where {q a } is a coordinate system in M./G. The kinetic part for the gauge-invariant per- 
turbations is [] 

l±Q a p±r±^ (el) . (3.22) 

Thus the form of ±Q a /3, which is necessary and sufficient to obtain the explicit form of the 
matrix ( [3.11| ), can be read from the kinetic part of the reduced Lagrangian. 

When the dot ' is identified with V F , L must derive the field equation ( |3.9|) with j_uf 
replaced by ±(p a since by definition the field equation ( |3.9|) is nothing else the evolution 
equation of the orthogonal perturbations. From this fact, the quadratic part of the reduced 
Lagrangian must be of the form 



^ a p±r±^-^±V a ^ a ±<P P (el) , (3.23) 



4 Note that terms linear in ^(j) a is automatically zero because of the background equation. 
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where ±V a p is defined by ( |3.15| ). The effective Hamiltonian ( |3.13|) can be obtained from the 



effective action by usual prescription. Thus, in a continuous limit, the reduced Lagrangian 
leads us to the effective quantum field theory considered in the previous subsection. 

After all, when we intend to investigate a state of the physical perturbations after a 
quantum tunneling, what we have to do is: (a) to fix a MPEP; (b) to give a coordinate 
system in a space of all gauge- invariant perturbations around the MPEP; (c) to obtain the 
reduced Lagrangian of the gauge-invariant perturbations by simply reducing the original 
Lagrangian to the space of all gauge-invariant perturbations; (d) to use the quantum field 
theory of the perturbations constructed from the reduced Lagrangian. From the effective 
action we can read off the form of the matrix (|3.11|) , which must be real and positive in the 
subspace generated by gauge-invariant perturbations in the Euclidean region. 

IV. SUMMARY AND DISCUSSION 

Throughout this paper, we have investigated the simultaneous differential equations ( |2.2[ ) 
with the algebra ( |2.5[ - |2.7| ) by using the WKB method. We have mentioned that, in a contin- 
uous limit, (|2.2j ) and (^5| - ^7| ) become generalizations of the Wheeler-DeWitt equation ( |1.8|) 
and the Dirac algebra ( |1.9| - |l.ll1) , respectively. We have defined an Euclidean wave function 



as that solution of the equations ( |2.2|) whose lowest-order part in the WKB expansion (|2.21|) 



is real in a region of the configuration space. We have called the region an Euclidean region. 
The lowest-order part in h has been expanded in the superspace around an instanton, by 
using a deviation equation of a vector field tangent to a congruence of instantons. The 
instanton around which we expand the wave function corresponds to a so-called most prob- 
able escape path (MPEP). Then we have shown that, when the expanded wave function 
is analytically-continued beyond the Euclidean region, the continued wave function can be 
understood in terms of a quantum-mechanical system, which will become a quantum field 
theory in a continuous limit. The corresponding state of physical perturbations around the 
analytically-continued classical path is the vacuum state determined by 'positive-frequency 
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mode functions' {_!_«£}. Each 'mode function' must be normalized with respect to the norm 
( p.4| ) and satisfy the evolution equations ( |3.9| - |3.10D , which correspond to linearized classi- 
cal equations of motion for perturbations orthogonal to the analytically-continued classical 
path. The positive-frequency mode functions must satisfy also the following boundary con- 
dition: the matrix Q a p defined by ( |3.11| ) must be real and positive definite in the subspace 
generated by {±Z?} when it is analytically-continued back to the Euclidean region. Note 
that the boundary condition of the mode functions is dependent of the matching condition 



of the WKB wave functions. In subsection [III A| , in the case that the weighted 'potential' V 



does not depend on (p n near the turning point with an enough accuracy, we have adopted 
such a matching condition that the WKB wave functions are analytic continuations of each 
other. If a future analysis for more general case gives a different matching condition of the 
WKB wave functions, then the mode functions must be matched with those in the Euclidean 
region in a different way and the corresponding matrix n a p must be real and positive definite 
in the subspace generated by {±Z?}. Anyway, the state of perturbations after a quantum 
tunneling is determined uniquely by the positive-frequency mode functions. Thus a quan- 
tum field theory is effective to investigate a state of physical perturbations after a quantum 
tunneling. We have shown that the effective Lagrangian describing the field theory is ob- 
tained by simply reducing the original Lagrangian to the subspace spanned by the physical 
perturbations around the MPEP. 

The result of this paper does not depend on the operator ordering since we have con- 
centrated only on ( [2.8| - 12.1lD and Q2.22|) , which are independent of the operator ordering. 



Moreover the result can be applied to a general background, which corresponds to a most 
probable escape path (MPEP). Thus, for a general MPEP, a quantum field theory is ef- 
fective to investigate a state of all physical perturbations after a quantum tunneling with 
gravity. Due to the result of this paper we can safely use the quantum field theory to obtain 
a CMB anisotropy, a spectrum of primordial gravitational waves, etc. in some inflationary 
scenarios. For example, in the one-bubble inflationary scenario, Tanaka and Sasaki [T^ 



calculated the reduced Lagrangian and quantized the gravitational perturbations. Although 
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they did not justify their treatment from the point of view of the Wheeler-DeWitt equation, 
the result of this paper do justify it without writing down the explicit form of the perturbed 
Wheeler-DeWitt equation. 

Finally we propose a possible definition of a two point function to strengthen the field 
theoretical interpretation. For a function / on Ai/G, define a one-parameter family of 
expectation values ( f )t along the analytically-continued classical path. 



</>* 



I(Und±^)J\det' ± g\\^f 



(4.1) 



t,V x =Vo 



where det' represents a determinant of the followed matrix restricted to the subspace gen- 
erated by {(d I d ±f n ) a } ■ It is evident that for the wave function ( |3.2| ) the expectation value 
of ±<f) a is zero: (±<f) a )t = 0. On the other hand the expectation value of {±<f) a , ±<fi } is not 
zero: 



({±r,±4> l3 })t-n[(n + ^)- 1 



a/3 



+ h 



0a 



(4.2) 



where the right hand side is calculated on the classical path. Note that means that 
the difference between the right and the left hand side is sufficiently small when the wave 
function is so peaked along the classical path that an orthogonal transformation of the 
hermite part of Q and the evaluation of the expectation value is approximately commutable. 
From (13.71) we can show that 



(Q + fi f ) 1 = ju + (_Lit, ±u) KG± u 



(4.3) 



and 



({±<p a , ±<P l3 })t ~ h [±u t {±u, ±u) K 1 G± u 



1 a/3 



+ h 



±u\±u, ±u) K 1 G ±u 



-i /3a 



(4.4) 



In subsection [LI [B| we have shown that the matrix Q: can be chosen so that (±Uk, ±Uk>)KG 
5kk'- For this choice of the two point function is reduced to 



(4.5) 
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The right hand side of (|4.5|) can be understood as a two point function of the quantum- 
mechanical system introduced in subsection as follows. The representation ( |3.16| ) leads 



the two point function 

(0\{J a ,±^}\0) = h^2(±nl a ± 4 + jL«Jx«f) , (4.6) 

k 

where the state |0) is the vacuum state of the quantum- mechanical system: 

Ofc|0) = (4.7) 



for v /c. The two point function ( |4.6|) is the right hand side of ( [4.5|) itself. Thus the two 
point function based on the wave function (|3.2| ) is equivalent to that based on the quantum- 
mechanical system. 
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APPENDIX A: DEVIATION EQUATION OF A VECTOR FIELD 

In this appendix we derive a deviation equation for a spacelike (and timelike) vector 
field in a pseudo-Riemannian manifold. By spacelike (timelike) we mean that its norm with 
respect to the metric is positive (negative). 

Let (M, g a p) be a pseudo-Riemannian manifold and N a be a vector field on M. It is 
assumed that g ,^N a N^ ^ in a region E (e M). Then introduce another vector field Z a 
in E such that 

[N,Z} a = . (Al) 

Locally this condition is a necessary and sufficient condition in order for N a and Z a to define 
independent coordinate variables by 
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- z ° ■ (A2) 

When we apply equations derived in this appendix to the system investigated in subsec- 
tion [II B| , N a corresponds to a gradient vector filed (|2.40|) of a real solution of ( [2.39|) . Then 



A corresponds to a coordinate ip n distinguishing physically different integral lines of N a . 
The condition (|A1|) can be written explicitly as 



Z a = Z p V p N a , (A3) 
where V is a covariant derivative compatible with g a p and 

X a = N p V p X a (A4) 
for X a G TM. Operating p to it, we obtain the following equation. 

Z a = zP\7pN a - R a pj3(T N p N a Z 13 , (A5) 

where R(X, Y)Z 

Without any modification, these equations are not effective enough to investigate the 
dependence of the wave function on the perturbation around the instanton. It is because 
Z a may have a component in the direction of N a . We want to modify these equations so 
that Z a appears as a vector projected to a hypersurface orthogonal to N a in the modified 
equations. To make the modification possible, we define a so-called Fermi transported basis. 
First a Fermi derivative VV along N a is defined by 

V F X a = X a + ±g{X, N)N a - ±g(X, N)N (A6) 

for X a E TM, where N 2 = g(N, N) Q. Note that N 2 is positive (negative) when N a is 
spacelike (timelike). The Fermi derivative has the following well known properties. 

1. If N a is tangent to a congruence of geodesies, then VfA° = X a for X a e TM. 
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2. V F± X a = j_( ± X)' a for X a G TM, where ± X a = X a - N a g(N,X)/N 2 . 



3. W F (N a /^/\N 2 \) = 0. 

4. If X a and F a (e TM) satisfy V F X Q = V F y a = 0, then N a d a g(X, Y) = 0. 

The first and the second properties suggest that the Fermi derivative is effective for our 
purpose. From the third and the forth properties, along an integral curve of N a the Fermi 
transported basis can be defined by 

e« = N°/J\W\ , 
V F ef = , (A7) 

and so that g(eo, e*) = and g(ei, Cj) = rjij, where ijij is a constant regular matrix. By using 
this basis, we can expand Z a as 

Z a = Z°e^ + Z i ef , (A8) 

where Z° = g(Z,e )/g(e ,e ) and Z* = (tT 1 )'* g(Z,e 3 ). Substituting flA8|) into flX3]) and 
we can show by explicit calculation that 

^-Z i = Z 3 V,N l , 
or 

d 2 ■ 1 • • 1 • 

Z l = Z^VjN* - E! ja N p N a Z J - —N l Z J r] jk N k - —N'Z^VjN 2 , (A9) 



q t 2~ 'J- P3<y- N 2 J N 2 

where 



V,-(iV 2 ) = e^iV 2 ) . (A10) 



The result can be rewritten covariant form: 
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where 



V 2 F± Z a = ± ZP ± (VpN) a - ± R a pPa NON° ± ZP - l^ g ( ± Z,N) ± N a 

-j^9(±Z,VN 2 ) ± N a , (All) 



N a N 
N a N 

(y pN y = VpN «-__L VpN i > 
N a N 1 

^R a pf3(T = -^"p/Jct J^^pPo- ■ (A 12) 
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